c3 



ON LIMIT CYCLES APPEARING BY POLYNOMIAL 
PERTURBATION OF DARBOUXIAN INTEGRABLE 

SYSTEMS 



DMITRY NOVIKOV* 



o 
o 

<N 

5-1 , 

Abstract. We prove an existential finiteness result for integrals of ra- 
tional 1-forms over the level curves of Darbouxian integrals. 

(N 

^ 1. Limit cycles born by perturbation of integrable systems 

1.1. Poincare Pontryagin integral. Limit cycles (isolated periodic tra- 
. jectories) of polynomial planar vector fields can be produced by perturbing 

integrable systems which have nested continuous families of non-isolated 
periodic trajectories. The number and position of limit cycles born in such 
perturbations is determined by the number and position of isolated zeros 
of the Poincare-Pontryagin integral of the dissipation form along the closed 
periodic orbits of the non-perturbed integrable vector field. The infinitesi- 
C^- \ mal Hilbert problem is to place an upper bound for the number of isolated 

zeros of this integral. This problem was repeatedly formulated as a relaxed 
form of the Hilbert 16th problem by V. Arnold [1]. 

Instead of polynomial vector fields, it is more convenient to deal with 
■ (singular) foliations of the real plane R 2 by solutions of rational Pfaffian 

! equations 9 = 0, where 9 is a 1-form on R 2 with rational coefficients (since 

only the distribution of null spaces of the form 9 makes geometric sense, one 
can always replace the rational form by a polynomial one) . The foliation T is 
(Darbouxian) integrable, if the form is closed, d9 = 0. An integrable foliation 
5_i ■ always has a "multivalued" first integral of the form f(x,y) = expr(x,y) • 

Pj{x, y) Xj , where r(x, y) is a rational function and pj polynomials in x, y, 
which are involved in the in general non-rational powers Xj. A particular 
case of the integrable foliations consists of Hamiltonian foliations defined by 
the exact polynomial 1-form 9 = dh, h € M[x,y]. 

If L C M 2 is a compact smooth leaf (oval) of an integrable foliation 9 = 
belonging to the level curve {/ = a} C M 2 , then all nearby leaves of this 
foliation are also closed and form a continuous family of ovals belonging to 
the level curves {/ = t}, t € (R 1 , a). The ovals from this family, denoted by 
Lt, are uniquely parameterized by the values of the real variable t varying 
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in an open interval a < t < [3. Their union is an annulus A*,/3 bounded 
by the Hausdorff limits L a = lim t _ >a , + o L t and Lg = lim t _ >/ g_o. The limits 
themselves may be not leaves of the foliation, but rather separatrix polygons 
(unions of several leaves and one or more singular points of the foliation J-). 
The polar locus of the form 6 is usually a separatrix polygon. 

Consider an one-parameter family of planar real foliations T e defined by 
the Pfafnan equations 

9 + euj = 0, Poles(w) C Poles(0), (1.1) 

with a closed rational 1-form 8 and arbitrary rational 1-form u. The as- 
sumption on the poles guarantees that the perturbation (1.1) does not create 
additional singularities in the annulus A a ^ (such u will be called admissible 
for 8). The Poincare-Pontryagin integral 1 associated with this family, is the 
function 

t ^ I{t) = j> w, L t C {/ = t} an oval, t £ (a, (3). (1.2) 

The following well-known Poincare-Pontryagin criterion establishes con- 
nections between the integral (1.2) and limit cycles of the system (1.1): If 
the function I{t) has an isolated zero of multiplicity fj, at an interior point 
a £ (a, f3), then the perturbed foliation T e for all sufficiently small values 
of e and 5 > has no more than fx limit cycles in the thin annulus between 
the ovals L a ^s and L a+ s- The assertion is reasonably sharp: under certain 
additional assumptions of nondegeneracy one can produce exactly fj, limit 
cycles. 

Thus bounds for the number of isolated zeros for the integral (1.2) trans- 
late into bounds for the number of limit cycles for near-integrable systems. 2 
The following conjecture seems to be believed by most experts in the area. 

Conjecture 1. For any pair of rational 1 -forms (9,lo) on the real plane M 2 
of degrees n,m respectively, such that dO = and Poles(u;) C Poles(#), the 
number of isolated real zeros of the integral (1.2) is bounded by a constant 
N = N(n,m) depending only on n and m. 

Demonstration (or refutation) of this Conjecture constitutes the infini- 
tesimal Hilbert Sixteenth problem [1,5]. A separate question concerns com- 
putability (and practical computation) of the bound N(n,m). 

1.2. Analytic properties of the Poincare Pontryagin integral: Ham- 
iltonian versus Darbouxian cases. The function I(t) defined by the in- 
tegral (1.2) is obviously real analytic inside the real interval (a, f3), and hence 
the number of isolated zeros is necessarily finite on any compact subinterval 

lr The function I(t) is sometimes called Melnikov function or (misleadingly) the Abelian 
integral, see below. 

2 However, the above Poincare-Pontryagin criterion does not apply to limit cycles born 
from separatrix polygons L a , Lp which bound the annuli of periodic orbits of the integrable 
systems. This bifurcation requires much deeper analysis that will not be discussed here. 
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of this interval. Yet this circumstance neither implies any explicit upper 
bound, nor rules out the accumulation of infinitely many isolated roots to 
the boundary points a, (3 of the interval. 

Most studies of the function I(t) are based on its analytic continuation 
from the real interval (a, (3) to the complex domain. Already at this step 
drastic differences between Hamiltonian and general integrable cases occur. 

1.2.1. Abelian integrals. If the initial foliation T = {9 = 0} is Hamiltonian, 
i.e., if 6 = dh, h G M.[x, y], then oj must also be polynomial and the function 
I(t) is an Abelian integral, integral of a polynomial 1-form over a cycle on 
algebraic curve. Such integrals possess the following properties. 

(i) The integral I(t) admits analytic continuation on the entire Riemann 
sphere CP 1 = CU{oo} as a multivalued function ramified over finitely many 
points. For a generic Hamiltonian h these points are exactly the (complex) 
critical values of h. 

(ii) I(t) is (a linear combination of coordinates of) a solution of a Fuchsian 
linear system on CP 1 . 

(iii) Near each singular point to of the Fuchsian system, I(t) admits a local 
representation as the finite sum I(t) = Y2j k a jk(t) {t — to) Xj ln k (t — to) with 
the coefficients ajk holomorphic at the point to and finitely many rational 
exponents Xj G Q. 

The last property alone suffices to guarantee that Abelian integrals are 
non- oscillating: any such integral has only finitely many real isolated zeros 
on any interval. This can be proved by application of the "derivation-division 
algorithm" [15]. Moreover, the Fewnomial theory [3] which can be regarded 
as a multidimensional generalization of the above algorithm, implies the 
following uniform version of the non-oscillation property, achieved in [4, 17]. 

Theorem 1 (A. Khovanskii-A. Varchenko, 1984). The total number of real 
isolated zeros of Abelian integrals over the level curves of a polynomial h 
is uniformly bounded over all polynomial 1-forms lo and all Hamiltonians h 
sufficiently close to any given combination (ho,Uo). 

By the standard compactness and homogeneity arguments, Theorem 1 
implies that for any combination of n, m £ N the number of real isolated 
zeros of Abelian integrals of forms of degree ^ m over Hamiltonians of degree 
^ n is bounded by a constant No(n, m) depending only on n, m. However, 
this constant is absolutely existential (non-constructive). Other, completely 
different tools allow to place explicit upper bounds on the number of zeros 
of Abelian integrals which are at least ^-distant from the critical values of 
the Hamiltonian [9, 10] (the bounds depend on S > 0). 

One can expect similar behavior of the Poincare-Pontryagin integrals in 
the case when 6 is exact albeit non-polynomial anymore, 6 = df, with 
f(x,y) = a rational function, though the subject was not explored 

systematically. 
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1.2.2. Pseudoabelian integrals. Completely different is the picture if the (un- 
perturbed) closed form 9 is non-exact. In this article we consider only the 
simplest case when 6 has at most first order poles after extension on CP 2 
and hence a Darbouxian integral of the form 

f{x,y)=l[ Pj (x,y) x i, Xj e R+, j = 1,. . . ,m, (1.3) 

j 

ramified over the separatrix polygon S(8) = \Jj{pj = 0} C K 2 corresponding 
to the "critical" value / = 0. For lack of better name, we will refer to 
integrals of the form 

I(t) = I u, L t C{f = t}C R 2 , (1.4) 
JL t 

where u is a rational form without singularities on Lj, as pseudoabelian 
integrals. 

This particular case already contains all difficulties which distinguish the 
general integrable case from the Hamiltonian one. The closed level curves 
{/ = t] cease to be algebraic unless all ratios \j/\k are rational, more- 
over, their complexifications are generically dense in C 2 [6]. The topological 
methods of continuation of the corresponding integral to the complex do- 
main, which work so nicely in the Hamiltonian case, fail [12], thus the ana- 
lytical continuation of the corresponding integral I{t) should be achieved by 
completely different methods. Besides, I(t) is not known to satisfy any rea- 
sonable linear or nonlinear differential equation of finite order, which makes 
impossible applications of the methods from [9, 10]. For example, the local 
representation (1.5) already implies that I(t) is not a solution of a Fuchsian 
ordinary differential equation of finite order. 

Even the most basic among the properties of the Abelian integrals, the 
local analytic representation, fails for pseudoabelian integrals. More pre- 
cisely, it can be shown that if the separatrix curve S(6) has only normal 
crossings and the collection of Darboux exponents [Ai : • • • : A n ] satisfies 
certain generic arithmetic properties, then the pseudoabelian integral I(t) 
can be locally near t = represented as a composition, 

I(t) = F(t,t^,--- ,t"»), lil ,..., fln €R+, (1.5) 

with a real analytic function F of n + 1 variables. This observation, due 
to C. Moura [8], suffices to prove non-oscillatory behavior of pseudoabelian 
integrals, but under the above arithmetical restrictions on the Darboux ex- 
ponents. Such form of arithmetics-conditioned finiteness is not unknown in 
the general theory of o- minimal structures, see [16]. Clearly, if this arith- 
metic restriction is indeed necessary to rule out non-oscillation of an individ- 
ual pseudoabelian integral I(t), there would be almost no hope for uniform 
bounds for the number of zeros of pseudoabelian integrals. 
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1.3. Principal result. Our main result is the uniform non-oscillation of 
pseudoabelian integrals associated with a generic Darbouxian integrable fo- 
liation near the separatrix polygon S(9). To formulate it accurately, we 
introduce Darbouxian classes of integrable systems. Let mi, . . . ,m n be a 
tuple of natural numbers, mj G N. 

Definition 1. The Darbouxian class 3> = f(mi, . . . , m n ) is the class of real 
integrable foliations on RP 2 defined by the Pfaffian equations 9 = 0, where 
6 is a logarithmic form with constraints on the degrees as follows 3 , 

9 = jr^Xj—, Pj G K[x,y], degpj < mj, Xj > 0. (1.6) 
j i Pj 

The polynomials are assumed to be irreducible and different, hence the 
polar locus S{9) is the union of n algebraic curves {pj = 0} C M 2 . 

Foliations from a fixed Darbouxian class are parameterized by points from 
an open subspace in the Euclidean space with the coordinates being the ex- 
ponents Ai, . . . , X n > and the respective polynomials pj, identified with 
their coefficients. Using this identification, one can define open neighbor- 
hoods of a given closed form 9 = in its Darbouxian class Q) . We will 
consider only foliations which have continuous families of ovals accumulat- 
ing to the separatrix S(9) or its proper subset (as usual, in the Hausdorff 
sense) . 

If a real rational 1-form to on R 2 is admissible for the logarithmic form 9, 
i.e., if the polar locus of to belongs to the separatrix S(9), then the integral 
I(t) = § L over any non-singular real oval Lt, t ^ 0, is well defined. We claim 
that the number of real isolated zeros of such integrals is locally uniformly 
bounded. 

Theorem 2 (uniform non-oscillation of pseudoabelian integrals). Let the 
form 9q = Y^j=i defining the Darbouxian integrable foliation T = 

{6*o = 0} be from the class S>(m\, . . . , m n ), and assume that the real algebraic 
curves {pj = 0} are smooth and intersect transversally (in particular, there 
are no triple intersections). Let loq be a rational one-form admissible for 9q 
with poles of order at most itiq. 

Then the number of isolated zeroes of the pseudoabelian integral near t = 
0, as a function of (9,uj), is locally bounded at (6q,ujq) over all admissible 
pairs (9, to). 

In other words, for any admissible pair of rational 1-forms (^O)^o) as 
in the Theorem, there exist a finite number N = N(9q,luq) and e > 0, 
depending on the pair, such that the number of real isolated zeros of any 
pseudoabelian integral corresponding to (9,u) and to a bounded family of 



One can easily show (e.g., see [6]) that a closed rational form 8 with poles of order 
^ 1 on CP 2 has the form (1.6) with complex constants \j. 
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cycles L t , t € (0, e), is no greater than N provided that the pair (0,lo) is 
sufficiently close to (6q,ujo). 

After this result one may expect that the number of zeros of pseudoa- 
belian integrals would be uniformly bounded over all Darbouxian classes. 
The following Conjecture is a relaxed form of the Infinitesimal Hilbert 16th 
Problem (Conjecture 1) 

Conjecture 2. For any Darbouxian class f (mi,...,m n ) and any degree 
mo there exists a finite number N = A" (mo, mi, . . . , m n ) such that any pseu- 
doabelian integral of a rational 1-form uo of degree ^ rriQ along ovals of the 
corresponding Darbouxian integral (1.3) has no more than N real isolated 
zeros. 

Again in contrast with the Hamiltonian case, the Conjecture does not fol- 
low from Theorem 2, since the compactness arguments fail for pseudoabelian 
integrals. Indeed, the assumption on the separatrix is an open condition in 
@i which in the limit may degenerate into non-transversal singularities at the 
corners. Besides, the Darbouxian classes themselves are non-compact and 
in the limit contain closed forms with poles of higher orders and "vanishing 
separatrices" as Xj — > + . All these scenarios require additional study. 

1.4. Proofs: Mellin transform and generalized Petrov operators. 

The most fundamental tool used to prove nonoscillation-type results, is 
the classical derivation-division algorithm [15] which goes as far back as 
to Descartes. In the modern language it consists in constructing a differen- 
tial operator T> (in general, with variable coefficients) such that the result 
of the derivation VI(t) is a function without real isolated zeros (e.g., iden- 
tical zero). By the Rolle-Descartes theorem, the number of zeros of zeros 
of I is then bounded by a constant explicitly expressible in terms of T> (in 
the simplest case of operators with constant coefficients, by the order of 
V, provided that all characteristic roots are real). This observation was 
generalized in several directions, among them for the vector-valued [3] and 
complex analytic [11] functions. 

However, this method encounters serious difficulties when applied to para- 
metric families of functions, especially the families which admit asymptotic 
(rather than convergent) representation as in (1.5) as t ^ + , with the 
exponents fij depending on parameters and eventually coinciding between 
themselves. The operator V, mentioned above, turns out to be very sensitive 
to these parameters. 

Our proof is based on introducing a new class of real pseutiodifferential 
operators V which decrease in a controlled way the number of isolated zeros. 
More precisely, for each such operator V and any family F of functions ana- 
lytic on the universal cover of the closed unit disk and analytically depending 
on parameter there exists a finite number ("Rolle index") A = N(V, F) such 
that for any u = u(t) € F on (0, 1) the numbers of real isolated zeros of u 
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and Vu are related by the "Rolle inequality" 

#{u = 0} ^ #{Vu = 0} + N. (1.7) 

These operators, depending on auxiliary parameters, in a sense interpolate 
between the "usual" differential operators, the Petrov difference operators 
that were used implicitly in [13] and explicitly in [14], and the operator of 
"taking imaginary part" which played the key role in [10]. 

The operator V is defined via the Mellin transform which associates with 
each function u(t) (say, defined on the interval [0, 1]), the function v(s) of a 
complex argument s by the formula 

v = Jtu, v(s) = [ t s ~ l u{t)dt. (1.8) 
Jo 

We show (Theorem 6 below) that the Mellin transform ^1 of a pseu- 
doabelian integral I(t) is a function holomorphic in the right half-plane 
Res > sq, which admits analytic continuation as a meromorphic function in 
the entire plane s € C with poles located on finitely many real arithmetic 
progressions on the shifted negative semiaxis M_. The distance between 
subsequent poles may be very small, depending on the arithmetic of the 
exponents Xj, which results in the divergence of the representation (1.5). 
The operator P is constructed as 

P = Jt~ 1 KJZ, K:v(s)^ K(s)v{s), (1.9) 

where K(s) is an entire function with zeros on these progressions. Then the 
product K ■ <s^I becomes an entire function whose inverse Mellin transform 
is identically zero. Note that a differential operator with constant co- 
efficients (the class sufficient for many applications of the derivation-division 
algorithm) can be represented in this way with a polynomial kernel K(s): 
this explains the connections with the traditional method. 

Flexibility of the construction of the operator V (location of zeros of K(s)) 
allows to choose V = Vq in a way analytically depending on the form £ 2$ 
such that V$I = for any pseudoabelian integral / = Ig(t) associated with 
the integrable foliation 6 = 0. 

The cornerstone of the proof is the demonstration of the Rolle inequality 
(1.7) for the operator V. We prove that the corresponding constant N = 
N(Ve) depends in a uniformly bounded way on the logarithmic 1-form 9 € 
Ql. This is sufficient for the proof of Theorem 2. 

1.5. Acknowledgments. I grateful to S. Yakovenko for drawing my atten- 
tion to this problem and numerous useful discussion. I am also grateful to 
C. Moura whose unpublished preprint [8] was the starting point of all my 
investigations. 
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2. Analytic representation of pseudoabelian integrals 

2.1. Parametrization of admissible pairs of forms. We define the com- 
plexification £Fc = ^c(mi,...,m„) of ^ as the space of logarithmic forms 



where, as above, pj are assumed to be coprime and irreducible. 

Denote by the subset f C f C £Fc of all 9 E ^ with smooth real 
algebraic curves {pj = 0} intersecting transversally. 

The pairs (9,u>) with 9 € and real one-form a; admissible for 9 with 
poles of order less than mo form a set = J (mo, mi, ...,m n ), which has an 
evident topology of a real analytic vector bundle with base . 

The set is a real analytic subvariety of the set B§Ci defined as the set 
of pairs of complex one- forms (9,uj), with 9 € 3>c an d u> admissible for 9 
with poles of order less than mo- The set has an evident structure of a 
complex analytic vector bundle over complex analytic manifold S>c- 

2.2. Local linearization of logarithmic foliations. We start with a well- 
known fact whose demonstration is only briefly sketched for completeness 
of exposition. 

Lemma 1. A closed meromorphic 1-form 9 with first order poles on two 
transversally intersecting analytic curves in C 2 can be analytically linearized 
near the point of their intersection. More exactly, there is a linearizing 
transformation <p mapping a neighborhood of a closed polydisk {\x\,\y\ ^ 
2} C C 2 to a neighborhood of the point of the intersection, such that the pull- 
back <p*9 is the logarithmic one-form A^r + A 4 "^ defined in this neighborhood. 
This linearizing transformation depends analytically on the complex form 



Proof. The two curves can be simultaneously rectified to become the coordi- 
nate axes. The local fundamental group of the complement is commutative, 
and is generated by two loops around these axes: this means that there 
exist two (complex) numbers A, fi (the residues of 9) such that the difference 
9 — (A^r + A*^) is exact, the differential of a meromorphic function. Because 
of the constraints on the order of the poles, this difference is the differential 
of a holomorphic function dg, holomorphically depending on 9. The inverse 
of the holomorphic transformation (x, y) i— > (x, y exp(/x _1 g)) brings the form 
9 to the required "linear" normal form in some neighborhood of the origin, 
which can be subsequently expanded to the polydisk by a linear dilatation 
in x, y. □ 




n 




9£@c- 
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Figure 1. Splitting of a cycle into pieces. 



Here is the linearization transformation written explicitly in terms of the 
first integral. Let 

n 

f = C Pl (x, y) x P 2(x, yY II Pj ( x < v) Xj (2- 1 ) 

i=3 

be a first integral of 9, C ^ 0, and assume that the transversally intersecting 
curves are {p\ = 0} and {p 2 = 0}. The inverse of the mapping 

n 

(x,y) ->■ (x,y), where x = Pl (x,y), y = p 2 Pj(x, y) Aj/M 

3=3 

is the required diffeomorphism 0. 

Lemma 2. For any C bigger than some Co 3> 1 the dilatation can be chosen 
in such a way that the first integral (2.1) becomes a monomial x x y^ in the 
new coordinates. 

Indeed, after linearization as in Lemma 1 the first integral becomes const x^y 11 , 
and if const > 1, then an additional dilatation brings the first integral to 
the required form. □ 

2.3. Representation for pseudoabelian integral. 

2.3.1. Local computation near saddles. We compute I{t) by evaluating sep- 
arately integrals along pieces Stj of Lt lying near saddles {pi = pj = 0} and 
along the pieces Sj of L t lying near smooth pieces of {pj = 0}. 

Integrals of monomial 1-forms along the piece 5ij of the integral curve 
L t = {/ = t} lying between two cross-sections, {x = 1} and {y = 1}, can be 
explicitly computed. Indeed, in the linearized coordinates (x, y) the curve 
5ij is, by Lemma 2, given by the explicit formula 



Sij = {y = ^x- A H t £ (0, 1), x £ (t 1 ^, 1), (2.2) 
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f 1 x p - 1 y q dx = t q/ » 



eventually after replacing / by Cf, with C> 1. Therefore the integration 
of uj = x p ~ 1 y q dx, (p, q) 6Z 2 , is simple: 

p-qX/H 1 fqjjx _ f p/\ 

— — rr- =^ n , , (2-3) 

P ~ qX/P t i/x p/X - q/fi 

if qj ' ji ^ p/X. For the resonant case q/p = p/X the same computation yields 
the answer 

i gt w hich can be otherwise obtained by direct passage 
to limit. Computation of f x x p y q ~ 1 dy is similar. 

Denote by £ pq x^(t) the following elementary function: 

' t p / x — t q ^ 

if^,/A# " / "• (2.4) 
fp/ x log t otherwise, 

One can verify that £ pq xn(t) depends analytically on A, p for any t > 0. 

2.3.2. Analytic representation of pseudoabelian integral. 

Lemma 3. After some rescaling of t, the pseudoabelian integral I(t) corre- 
sponding to (6,oj) sufficiently close to (9o,ujq) in£$ admits the representation 
near the point t = 

J(t)= + J>rA* r/ \ ( 2 -5) 

p,q,X,fi r,X 

where p, q, r run over a translated octant Z,^_, p,q,r € — tuq + N, and the real 
indices A, p run independently over the set A = {Ai, . . . , A n } C R™ of the 
Darboux exponents of the closed 1-form 6. 

Moreover, the coefficients a pq \^,b r \ depend analytically on the forms 6,u>, 
and can be continued holomorphically to some neighborhood U of (6q,ujo) G 
£$c, an d continuously to its closure U <s S§c- The coefficients a,A M ,6#A 
decrease at least exponentially: there exist C, p > 2 such that the following 
bounds hold in U 

\apqx»\^Cp-( p+q \\b rX \^Cp- r . (2.6) 

Proof. By assumptions, the Hausdorff limit Lq = lim t _ >0 - L t C S{9) of the 
real ovals is a bounded curvilinear separatrix polygon, whose edges and 
vertices depend analytically on 6 £ £F in a sufficiently small neighborhood 
of #o by the Implicit Function Theorem. 

Near each vertex one can linearize the foliation using Corollary 1 in such 
a way that the first integral is a monomial in the new coordinates, / = x A y M , 
eventually after a rescaling of /. 

The linearizing transformation brings the rational form uj into a mero- 
morphic form in the bidisk {\x\, \y\ ^ 2} with poles of order less than mo on 
the coordinate axes. Such form can be represented as a series 

^ a pq x p ~ 1 y q dx + (3 pq x p y q ^ 1 dy 

p,q>-m 
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converging in the bidisk {|x|, |y| ^ 2}. This convergence implies the expo- 
nentially decreasing upper bounds (2.6) on the Taylor coefficients a pq , (3 pq 
of the form lo in the linearizing coordinates x, y. Since the linearization 
depends analytically on 9 in a closure U <s of a sufficiently small neigh- 
borhood U of #o G by Lemma 1, the coefficients a pq ,f3pq are analytic 
in U and the upper bounds are uniform in this neighborhood. Integrating 
the series termwise along 5ij(t) = {/ = t} n {\x\, \y\ ^ 2} as in (2.3), and 
summing up over all vertices of L corresponding to the pair (A, //) (their 
number is finite by Bezout theorem) we obtain the first sum in (2.5) and 
the exponential upper bounds (2.6) for a pq \^, as finite linear combinations 
of oi pq , (3 pq . 

The local cross-sections near the corners given in linearized coordinates 
by {x = 1} and {y = 1} become analytic curves transversal to S(8) in the 
initial coordinates on the plane. The integral of the form lo along arcs of 
leaves L t between the two cross-sections near two endpoints of the j'th edge 
(on the curve {pj = 0}), depends analytically on the arc, as the latter re- 
mains nonsingular. In the chart t obtained by restriction of the Darbouxian 
integral, this means that the corresponding contribution to the integral is 
an analytic function of t 1 ^ 3 , where Xj is the residue (Darboux exponent) 
corresponding to the edge. Again, after a rescaling one can assume that this 
analytic function converges in a small nei ghborhood of {(t 1 ^) <j 2 }, uni- 
formly over a small neighborhood of (9q,ujo) G £%c, thus giving the second 
inequality in (2.6). □ 

3. The class Z and its Mellin transform 

We introduce a class Z of functions of one variable as sums of series of 
type (2.5) with coefficients satisfying (2.6), and define the notion of analytic 
3-family in Section 4 using the statement of Lemma 3 as its definition. We 
prove in Section 5 that the number of zeros on [0, 1] is uniformly bounded 
for any analytic ^-family. 

3.1. Definition of the class Z- 

Definition 2. We define ^-series as a formal series a of the form 

^ a pq\^pq\pi{t) + ^ b r \t T ^ X , UpqXfi, K\ G C, (3.1) 

p,q,\,V r,A 

where A, \i G A = {Ai, A n } C M+ - a finite set, p,q,r <G — m a + N, and 

\a pq \»\ < Cp- p - q , \b rX \ < Cp~ r for some C,p > 2. (3.2) 
The set A is called the spectrum of a. 

We denote further by M = M(a) some number smaller than all exponents 
in (3.1), e.g.: 

M = mm(-m a \i). (3.3) 

i 

The class Z is the class of functions f(t) in one variable t G [0, 1], called 
3-functions, representable as a sum of some ^-series. 
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One should emphasize that the ^-series corresponding to a function / G (5 
is not unique: indeed, even the function (. pq \^{t) can be written as a sum of 
two monomials. Also, we do not know any effective way to check if a given 
function belongs to Z- However, according to Lemma 3, the pseudoabelian 
integrals belong to 3- 

3.2. Analytic continuation of 3-functions. We prove here that ^-functions 
can be analytically continued to the universal cover of the punctured unit 
disk. 

We start with an elementary inequality. 



Lemma 4. For any t, a, (3 G C we have 

t a -t? 



a- f3 



^ 1^1 1 log t| 



e az _ e /3z 




pa 




z / e wz dw 






J/3 



for some 7 G [a, (3] C C. 

Indeed, take z = log t in the following inequality 

^ \a — [3\ ■ |e 7Z | • \z\ for some 7 G [(3, a}. □ 

Let C* denote the universal cover of C* = C \ {0} and let D° C C* be 
the universal cover of the punctured closed unit disk. 

Lemma 5. For any A > the series (3.1) converges uniformly after mul- 
tiplication by t~~ M in a neighborhood in C* of the sector Sa = {\t\ ^ 
1, 1 argt| ^ A} C D° , and, in particular, on [0, 1]. 

Lemma 4 implies that for some constant Ca and for all p, q, A, in (3.1) we 
have \t" M e pq x^(t)\ ^ C A 2 (p+q) as long as \t\ l / x , \t\ 1 /' x < 2. These inequalities 
define a neighborhood of Sa in which the series (3.1) converges uniformly 
by (3.2). □ 

Corollary 1. ^-functions are real analytic functions on (0, 1] admitting an- 
alytic continuation to the universal cover of the punctured closed unit disc. 

Remark 1. If the ratios of all consecutive Darboux exponents Aj/Aj, i,j = 1, . . . ,n 
are either rational, or "nice" irrational (not admitting too close approximations 
by rational numbers so that the differences p\i — q\j decrease no faster than ex- 
ponentially with p,q — > +00), then the series (2.5) can be re-expandcd in powers 
of t 1 /^, and this re-expansion converges. Then the integral I(t) becomes repre- 
sentable under the form (1.5). Application of the standard fewnomial technique 
allowed C. Moura to prove that for such collections of the Darboux exponents the 
pseudoabelian integral I(t) is non-oscillating, moreover, that the number of isolated 
zeros is locally uniformly bounded over u>, provided that the Darbouxian exponents 
are fixed [8]. However, variation of the exponents clearly destroys the arithmetic 
conditions. 
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Remark 2. The transcendental binomials £ pq \ IJ , which in the limit converge to log- 
arithms, first appeared at least as early as in 1951 in the work by E. Leontovich [7] 
on bifurcations of separatrix loop. Since then they regularly re-appear in different 
disguises, most recently as the so called Ecalle-Roussarie compensator. 

This ubiquity can be probably explained by the fact that the function £ e {t) = 
t(t 6 — l)/e is in a sense the only possible Pfaffian deformation of the "logarithm" 
£o{t) = —tlogt. Indeed, the graph of the latter function is an integral curve of 
the linear equation £o = 0, where £o = tdy + (y — t) dt is a Pfaffian form. This 
linear Pfaffian equation has a singularity of the Poincare type (resonant node) at 
the origin. Any (nonlinear) deformation £ e of the form £o is necessarily analytically 
linearizable by the Poincare theorem, the linearizing chart depending analytically 
on the parameter e. Among integral curves of the perturbed linear system one is 
the graph of £ e . 

3.3. Mellin transform of %. The (one-sided) Mellin transform of a func- 
tion u(t) G Lj oc ((0, 1]) denned on the interval (0, 1] is the function v = 
of complex variable s defined by the integral transform (1.8): 

v (s) = Jtu{s) = f t?- 1 u(t)dt. 
Jo 

If the function u grows moderately as t — > + , \u(t)\ = o(t~ M ) for some 
finite M, then the integral in (1.8) converges uniformly on compact subsets 
of the half-plane {Res > M}. In particular, Mellin transforms are defined 
for all ^-functions. 

The formal Mellin transform of a ^-series is defined by its action on the 
basic functions 

•*(*' /A ) = JTrA' ■""»■*< ~ (. + P /A)'(. + «/<.) 

and extended by linearity, so the result is a formal sum of these rational 
functions. We prove that this formal sum converges to a meromorphic func- 
tion which is an analytic continuation of the Mellin transform of the sum of 
the 3-series. 

Lemma 6. The Mellin transform g = of a function f G 2 given as a 
sum of the series (3.1) is the sum of the series: 

^E ( ,^+^) + £lT7A- (3 ' 4) 

p,q,\,H r,X 

with the same p, q, r, A, ji and a pq x^, b r \ as in Definition 2. 

In particular, g(s) extends analytically as a meromorphic function on 
the entire complex plane C, the poles of g(s) are of at most second order, 
are all real and lie in the union of n arithmetic progressions of the form 

uf = A r Vo-N). 

Proof. The Mellin transforms of each elementary binomial £ pq x^ and of t r ^ x 
are equal to its formal Mellin transform. Termwise integration of right-hand 
side of (3.1) multiplied by i s_1 is possible due to the uniform convergence on 
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[0, 1] of this product for Res > — M by Corollary 5. Therefore we get (3.4) 
in {Res > — M}. This double series converges to a meromorphic function 
on C given by the same formula. □ 

Let 7m be the boundary of the semistrip {z\ Re z ^ — M+ 1, j Im z\ ^ — 1} 
oriented counterclockwise, with M defined in (3.3). 

For any function bounded on 7m we define an integral transform Cm as: 

C M g{t) = ^- j^t- s g{s)ds. (3.5) 

For a meromorphic in C function g bounded on {| Ims| = 1} with poles 
on the real line and bounded from the right we define Cg as liniM^-oo C-M9- 

Theorem 3. The Mellin transform of any function f G Z is bounded 
on 7m, where M = M(o) for some Z-series representing f . 

The restriction of the integral transform C to JtZ is the inverse operator 
to the Mellin transform ^ . 

For x £ R and f G Z the transform C of e~ t>cs \M f '(s) is defined and is 
equal to f{e l>( i): 

C(e-^ s ^f(s)) = f(e^t). (3.6) 

Abusing the language we call C the inverse Mellin transform, though its 
domain of definition is bigger than ^Z- 

Proof. Evidently |s + p/X\ ^ 1 on 7, so the convergence of (3.4) is uniform 
and \g\ ^ J2 P , q ,\,v \ a P^A + E r ,A \ h r\\ < 00 on 7. 

The kernel t~ s decreases exponentially on 7 for \t\ < 1, so we can integrate 
t s g(s) termwise on 7. Calculating the residues, we see that 

t ^ 1 C ( - : , „ s , : ; r ) = ^pqX^i-, 



J + r/XJ ' \{s+p/\){s + q/fi), 

so the Mellin transform and the inverse Mellin transform translate (3.1) to 
(3.4) and vice versa, i.e. are mutually inverse if restricted to Z and ^Z 
corr esp ondingly. 

Since e~ t>cs is bounded by e'^ on 7, one can compute the inverse Mellin 
transform (3.5) termwise, again due to the uniform convergence of the series 
for g(s) and the exponential decay of t~~ s . For £ pq Xfj, and t a the (3.6) can be 
checked immediately. 

□ 

3.4. Asymptotic series for ^-functions and quasianalyticity of Z- 

Let / G Z be a sum of some 3-series a. Any ^-series a can be expanded into 
a formal series by powers of t: 

f = E ta i^ 1 + c « 2 M ' a G (- m - + N ) A_1 ' ( 3 - 7 ) 

a 

where A is the spectrum of a and A -1 denotes the set of its reciprocals. We 
prove below that the result of the expansion does not depend on the choice 
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Figure 2. Splitting of the integration contour 7 into two. 



of a and represents the asymptotic series for /. Moreover, the class 3 turns 
out to be quasianalytic at 0: any two functions having the same asymptotic 
series coincide. 

Lemma 7. The coefficients c" 1 , c~ 2 are the Laurent coefficients of the Mellin 
transform g = f of f at the pole s = a. 

The series f is an asymptotic series for f, so does not depend on a. If 
f = then f = 0. Equivalently, the only entire function in is zero. 

Proof. The first claim follows because the Laurent coefficients of g are given 
by exactly the same formulae as c a ,d a : 

c q X = b w + X v/X-a/u = ^ apqX ^ 

<7,M p,q,\,ti 11 Hlh ^ p,q,Xn 

q/fi—a qlfji^zp/X =a q/fi=p/X=a 

both converging due to (3.2). 

Choose some < Ao ^ min A. Each segment [kXo, (k + l)Ao), k G Z, has 
at most n points common with AZ, so one can choose Mk from this segment 
at least Ao/2n-distant from AZ. 

This implies that \z + p/X\ ^ Xq/2u on Rez = Mk, so \g\ ^ C on the 
boundary 7^ of the semistrip {Re z ^ Mk, \ lmz\ ^ 1}, and the main point 
is that C is independent of k. Denote by 7^ the boundary of the rectangular 
{Mk ^ Rez ^ — M + 1, 1 lmz\ ^ 1} and let split integration along 7 in the 
definition of C into a sum of two integrals, one along 7^ and another along 
7j.j see Figure 3.4. 

Evaluation of the integral along 7I gives a finite partial sum of /: 

2^ / t- s g(s)ds = X) ( Ca + da log £ ) ' 

m J'Vk a<~M k 

Integral along 7^ can be estimated from above by 2Ct~ Mk (1 + 1/logt), 
which proves the second claim. 

If / = 0, then g has no poles at all, i.e. g is entire. Therefore the 
integral along 7I is zero for all k. Therefore \f(t)\ ^ 2Ct~ Mk (1 + 1/logi) 
for all t € (0,1), and for all Mk- Since Mk — » —00 as — > —00 and C 
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is independent of k, we conclude that f(t) = on (0,1), and therefore 
everywhere by analyticity. This means that g = f = as well. □ 

4. Analytic ^#3-families and operator 

Definition 3. Let V (s R n be some compact set, and let U C C n be a 
bounded neighborhood of V. We denote by DvcU the space of real analytic 
functions on V which can be extended holomorphically to U and continu- 
ously to the closure U of U, equipped with the norm \\u\\q := m ax^ g jj |it(f |). 

Lemma 8. The space DvcU is complete. □ 

Definition 4. Consider a family F = {ov} of ^-series, v € V <s W N , and 
denote by A(u) = {Xi(u), X n (u)} C M+ the spectrum of the 3-series cr(u). 

We say that F is an analytic family of ^-series on pair V C U for 
some bounded neighborhood U C of V in if the functions Aj(z/), 
a pq\(v)n(u){ 1 ')-, b r \(v)( u ) belong to DvcU, the analytic continuation of \i{y) 
does not vanish in U (i.e. \~ l {v) G Dycf/ as well), and 

\Kw*») WIId l < C7p _(P+,) , || 6rA M < ^P~ r (4-1) 

for some C = C(F),p = p(F) > 2. 

An analytic ^-family F = {/(£; ^)} C 3 is a family of sums of an analytic 
family of ^-series. In other words, for every f(t; v) one can choose a -3-series 
o~ v such that f(t; v) is a sum of o~ u , and {cr u } is an analytic family of ^-series. 

Analytic ^#-3-family is by definition the Mellin transform of an analytic 
-3-family. 

Lemma 9. For a given analytic family F of ^-series the number m a is 
uniformly bounded from above by some number mp < oo for all a = u v G F . 

Indeed, the sets Sk = \y € U\m Uu ^ k <G N} are analytic subsets of U 
and their union is U, so Sk = U for some k. □ 

Lemma 3 says that the family of pseudoabelian integrals parameterized 
by SB is locally an analytic ^-family. 

Let F = {a u } be an analytic family of ^-series on pair V C U, and consider 
the series 

a pqX{v)lx{y) l pqX{y)il{v) {t) + J2 h r\{u)t r/m leave U. (4.2) 
p,q,\(v),n(y) r,\(y) 

Lemma 10. For every A > there exists some e > such that this series 
converges uniformly on compact subsets of 

U £ = {| argi| ^ A,0 < \t\ < e 6 } x ({(ImAr 1 ) < e} n U) C C* x U, 

so its sum is a function holomorphic on U e . 

Proof In U £ we get |t A ~»| < exp (e [A + ||A _1 (^)||o] ), i.e. \£ pq x^(t)\ can 
grow exponentially fast. However, the fast decay of coefficients a pq \^,b r \ 
guarantees convergence if e is sufficiently small: the series (3.1) converges 
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uniformly on compact subsets of U e as soon as e [A+ || XT 1 {y) \\q] < In 2 for 
i = l, n. □ 

Remark 3. For functions from ^Z one can try to define analytic dependence 
on parameter as analytic dependence of f(t; v) on v on compact subsets of 
C\R. Analytic ^3-faniilies depend analytically on parameter in this sense, 
as shown above, but the opposite implication seems to be wrong. 

4.1. The operator E^. From Lemma 11 below follows that the rotation of 
argument preserves Z- for any f(t) G Z the function f(e l>t t), x £ M, is also 
in Z- From (3.6) it seems that the Mellin conjugate to the rotation is the 
operator of multiplication by e~ l>cs . However, the latter does not preserves 
^Z'- the function e~ ll<s ^ f (s) grows as | Ims| — > oo, while any function in 
^Z decreases as |Ims| -1 . 

In this section we define the Mellin conjugate £^ of the rotation on ^Z, 
and prove that application of operator £^ with x analytically depending on 
parameters preserves ^#3-families. 

Let start from the standard computation applied to / G Z- 

JZ(f(e i>e t))= f t s - 1 f{e i>c t)dt 
Jo 

= P t s - 1 f{e i>c t)dt+ [ t s - l f{e ll< t)dt 
Jo Jc* 

= e~ il<s C u s - l f{u)du + R{s) = e~ ii<s Jtf{s) + R(s), (4.3) 
Jo 

where = {e tz ,z G [— x, 0]} and R(s) = f c t s ~ 1 f(t)dt is an entire func- 
tion of exponential type. One can immediately see that R(s) is bounded in 
{|Ims| 1}. Therefore CR = 0, which explains (3.6). We take this as a 
definition of £^. In a sense E^ is an operator of projection of e~ l>cs ^Z to 
^Z along the space of entire functions. 

Definition 5. Let x G R. For any g G ^Z define E^g as the only function 
in ^Z such that the difference e~ li<s g — £n(g) is an entire function. 

Remark 4. One can show that £ >( g(s) is an analytic continuation from 
{Ims > 1} of the integral 

1 f e~ iitz g , 
-dz, 



2iri J 1 z - 
where 7 is as in (3.5). 

Lemma 11. Let F = {g(s;i')} be an analytic ^Z-jamily on -pair V C U, 
and let x = x(y) G QycU- Then 

(1) {£^t u )g(s; u)} is defined and is an analytic ^Z-family on pair V C U 
for some U C U. 
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(2) the operator C o o ^ maps analytic ^-families to analytic 2- 
families, 

(3) CSttJKf = {e ix t), and 

(4) the operator of rotation of argument preserves Z- 

Proof. The uniqueness in Definition 5 follows from Lemma 7: let g\ , gi G 
be two functions which differ from e~~ %7<s g by entire functions. Then 
their difference is entire and lies in J%%, which, by Lemma 7, implies that 
the difference is zero. 

Let us construct £ >c . This is simple for ( s+p/x ^ s+q/ ^ and ^p^: 

g— ixs ^—i?cx ^—ixs g— ixx 

= + , 

s — x s — X s — X 

and 

e -ixs e -ixx _j_ e -ixy e ~ixx _ e ~ixy / j j \ 

+ + + R, 



(s — x)(s — y) (s — x)(s — y) x — y \s — x 

(s— x)e~ ^^y+(x v\e~ i>cs -\-(v s)e~ ™ x 

where R = 2'- >- is \/ s L y){x „yf is entire in s. 

We extend this by linearity to any g € ^#3- Namely, for g given by (3.4) 
we define £^g(s) 

by the following formulae: 

<W = \ {^ /X + « M A„ (4.4) 

-ixr/X _ p-ixq/n 

Kx = e'-/\ x + £ — _ a rgA ,. (4.5) 

Take e so small that p = pe~ £ > 2. Let a neighborhood U C U of V 
be so small that |lm (x(i/)A i _1 (i/)) | < e for i = l,...,n and v € £/ '. Then 
| e -i*r/A| < e £r in U. 

For p,q,r > —m as in (3.4) we have max(p, q) < p + q + m, so for u G J7 
we have estimates 



2 



e -ixr/X _ e —ixq/n 



r/X - q/p 



by Lemma 4. 

The upper bounds (3.2) then imply that (4.4), (4.5) converge in D yc ^, 
and 

which proves the first claim of the Lemma. The second claim is a conse- 
quence of the first and of definition of ^3-family. The third claim follows 
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from (4.3) and (3.6), so the fourth follows from the third and the fact that 
Ex preserves J%Z- □ 

5. PETROV OPERATOR AND PROOF OF THEOREM 2 

We generalize the so-called Petrov operator, see [14], as follows: 

T«f{t) := i (f(te~n - f(ten) x G R (5.1) 
Lemma 12. The inverse Mellin transform of g(s) sin(xs) is V^f. For any 

v„f = j(~ x o (£„ - o ^/ 

This is a straightforward applications of (3.6) and of Lemma 11 corre- 
spondingly. □ 

The key idea of the proof is a generalization of the Petrov trick invented 
in [13], where it was used for x = it. 

Let f(t) be a ^-function. Denote by N(f) the number of isolated zeroes 
of function / on (0, 1) counted with multiplicities, by A^(/) the increment 
of arg/(t) along the arc {e"*^, <fi G [—1,1]} and let A^(/) be equal to the 
limit as e — > of the increment of arg/(t) along the arc {ee*^, G [—1, 1]}. 
Finiteness of all these numbers follows from Lemma 7. 

Lemma 13. For f G Z and real on (0, 1) the following inequality holds: 

N(f) < 1 + iV(7W) + + A° (/)). (5.2) 

Proof. Applying the argument principle to the sector {| arg z\ ^ x, e ^ |z| ^ 
1} on the universal cover of the closed unit disc we get 

2nN(f) < A_ + A + + A^/) + A° (/), (5.3) 

where A± are the increments of the argument of f{t) along the rays I± = 
{te ±l>( , t G [e, 1]}. Since f(t) is real on the real axis, it takes conjugate values 
on I±, so A + = A_ and Im/|/ ± = ±7^/. 

First, assume that V^f ^ 0. Evidently, any segment of I + where incre- 
ment of argument of f{t) is greater than tt contains at least one zero of 
Im/(t), so 

A + < vr (7V(Im f\ I+ + l) = nN{P„f) + vr, 
and we get the (5.2). 

If V x f = 0, then, by symmetry, f{te iie ) = f{te' i>c ), i.e. / = f(t n ^) with 
/ holomorphic in a punctured unit disc. Therefore (5.2) follows from the 
argument principle applied to / and a punctured unit disc. □ 

Remark 5. As x — > 0, P^f(t) tends to a multiple of f'(t), the contribu- 
tions A^(/) and A%(f) of arcs tend to zero, and we get the standard Rolle 
theorem as a limit case of Lemma 13. 

Let F = {f(t, u)} be a CJ-family on pair V C U, and let x = x(z/) g £Vct/- 
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Lemma 14. The variation A^^f is uniformly bounded from above by some 
A 1 = A 1 (F, x) < oo for all v G V. 

Proof. The functions Im. f{e i ^ <t> ;u) and Re f \e {l < v ^ '■ u) are analytic 
functions of (f> G [—1,1] analytically depending on parameter v G V <s 
M n . Therefore their number of isolated zeros counted with multiplicity is 
uniformly bounded from above by some number by Gabrielov's theorem [2]. 
The variation A 1 is at most 7r(this number +1). □ 

Lemma 15. The variation A^.^/ is uniformly bounded from above by some 
A = A°(F, x) < oo for all ueV. 

Proof. For a sum / £ 3 of a ^-series a the number A^.(/) is bounded 
from above by — m CT x by Lemma 7. Therefore the existence of the uniform 
upper bound follows from the Lemma 9. □ 

Corollary 2. For F, x as above and for any f G F 

7V(/) < 7V(P„/) + A(F, x), A(F, x) = A°(F, x) + A^F, x) + 1< oo. 

(5.4) 

5.1. Proof of Theorem 2. Starting from the analytic 3- family F of pseu- 
doabelian integrals parameterized by (6 : lj) G U C we construct the new 
family {P 7T x n ( I/ )f(», ^)}- It is again an analytic ^-family by Lemma 11. By 
Corollary 2, application of the operator P 7T \ n ( u ) reduces the number of iso- 
lated zeros at most by a finite number uniformly over parameters. 

On the other hand, the new family is simpler: its Mellin transform now 
has poles on {m/\k{v),k = 1, ...,n — l,m G Z>_ m/ }, i.e. on the union of 
only n — 1 arithmetic progressions. Indeed, Mellin transform of {P 7r \ n ( u )fu} 
has the same poles as sin(-7rA n (^)s)^#/, and the first factor has simple zeros 
at {m/\ n {u),m G Z>_ m/ }. 

Repeating this step n times we arrive at the analytic ^-family whose 
Mellin transform consists of entire functions. This means that this analytic 
^-family is identical zero by Lemma 7, so has no isolated zeros. □. 
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